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We examine effects of the Generalized Uncertainty Principle, predicted by various theories of 
quantum gravity to replace the Heisenberg's uncertainty principle near the Planck scale, on post 
inflation preheating in cosmology, and show that it can predict either an increase or a decrease in 
parametric resonance and a corresponding change in particle production. Possible implications are 
considered. 



I. INTRODUCTION 



It is generally believed that there was a so-called pre- 
heating phase near the end of inflation, in which energy 
was rapidly transferred from the inflaton to matter fields. 
This was followed by reheating, in which thermalization 
took place, and most of the standard model particles in 
our universe were produced [1-3]. During preheating, 
coherent oscillations of the inflaton field (f(t), around 
the minimum of its potential effectively contributed to 
a time-varying frequency term to the equations of mo- 
tion of matter fields x(*) coupled to it, thereby inducing 
instabilities by a well understood process known as para- 
metric resonance (PR) [4-7]. These in turn resulted in 
an explosive particle production [8-10] 1 (see also [1] for 
a recent review). Since the energy density of matter and 
radiation is exponentially small near the end of inflation, 
it has been argued that this transfer of energy from infla- 
ton to matter fields must be fast. It was recently shown 
that (in)homogeneous noise, such as those arising out 
of quantum fluctuations, indeed increased the instabil- 
ity band, and the rate of particle production [7] . Also it 
was shown in the past that these resonance effects were 
sensitive to non-linearities in the equations of motion for 
X, which in general were expected to increase the rate of 
particle production and result in an early termination of 
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1 It is also possible PR to result in a particle production which will 
not be explosive [11, 12]. We thank the anonymous referee for 
pointing this out to us. 



PR [5]. 

In this paper, we completely analytically study the 
effect of one such important non-linearity, that pre- 
dicted from the so-called Generalized Uncertainty Princi- 
ple (GUP). This in turn has been predicted from various 
approaches to Quantum Gravity, to replace the famil- 
iar Heisenberg's Uncertainty Principle near the Planck 
scale. We observe that, depending on the form of the 
GUP chosen, and initial conditions, an enhancement of 
particle production and an early termination of PR can 
indeed result. Our paper is organized as follows: in Sec- 
tion II we give a brief review of GUP. In Section III, we 
review relevant aspects of PR. Section IV considers ef- 
fects of general non-linear terms on PR and Section V 
the specific ones due to GUP. Section VI generalizes our 
analysis to an expanding universe. We conclude in Sec- 
tion VII. 



II. GENERALIZED UNCERTAINTY 
PRINCIPLE 



Various approaches to quantum gravity, including String 
Theory, Loop Quantum Gravity, as well as Black Hole 
Physics predict a minimum measurable length of the 
order of the Planck length, and a modification of the 
Heisenberg Uncertainty Principle to a so-called Gener- 
alized Uncertainty Principle, or GUP, near the Planck 
scale [13-19]. Such a form of GUP and the corresponding 
modified commutators between position and momentum 
coordinates was suggested earlier by two of the current 
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authors [20, 21] 



Ax t A Pl > -[l + l3((A P ) 2 +<p> 2 ) 

+ 2/3(Ap 4 2 +<K> 2 )] , » = 1,2,3 (1) 



[xi,pj] = ih(Sij + P&ijp 2 + 2j3piPj) 



(2) 



where p 2 = Ep 3 Pj> P = f3o/(M Pl c) 2 = P £ 2 Pl /2h 2 , 

3=1 

Mpi = Planck mass, and Mpic 2 = Planck energy w 
1.2 x 10 19 GeV. Note that the new (/3-dependent) term 
is quadratic in the Planck length (inverse Planck mass). 

Recently, some of us proposed another version of GUP 
consistent with the above as well as with Doubly Special 
Relativity (DSR) theories [22, 23]. Some phenomenolog- 
ical implications of the above GUP were also examined 
[22-25] This new version of GUP is of the form 



AxAp > - [1 - 2a < p > +4a 2 < p 2 >] 



> 




+ 4a z )Ap 2 +4:a 2 (p) 2 -2ay/{p 2 } 



[Xi,Pj\ 



ih Sij-a^pSij + + u 2 (p 2 Sij + SpiPj) 



(3) 



(4) 



where a = a /M P ic = a a £ P [/h. Note that in this case, 
the new (a-dependent) term is linear in the Planck length 
(inverse Planck mass). 

It is also noteworthy that although Eqs. (3) and (4) 
are not Lorentz covariant, they are DSR covariant [18]. 
We expect the results of our paper to have similar co- 
variance as well. In addition, since DSR transformations 
preserve not only the speed of light, but also the Planck 
momentum and the Planck length, it is not surprising 
that Eqs. (3) and (4) imply the following minimum mea- 
surable length and maximum measurable momentum 



Ax > (Ax), 
Ap < (Ap) r , 



a £pi 
Mpic 
a 



(5) 
(6) 



It is evident that the GUP induced terms in both (2) and 
(4) become important near the Planck scale. Also they 
both ensure, via the Jacobi identity, that [xi,Xj] = = 
\PhPj]- 

Next, defining 

Xi = xm , Pi = Poi (l + PpI) (7) 
and Xi = x Qi , p { = p oi (l - ap Q + 2a 2 pl) (8) 

3 

where p\ = ^ PojPoj and xoi , Poj satisfying the canoni- 

i=i 

cal commutation relations [a;oi,p<y] = *^ h is eas Y to 
show that equations (2) and (4) are satisfied to order /3 
and a 2 , respectively. Here, poi can be interpreted as the 
momentum at low energies (having the standard repre- 
sentation in position space, i.e. poi = —ihd/dxi), and pi 
as that at higher energies. 



Finally, employing Eqs. (7) and (8), it is easy to see that 
a non-relativistic Hamiltonian of the form 



H 



translates to 



7^+V(r) [f=( Xl ,x 3 ,x 3 )} (9) 



H = Ho + ^+OiP 2 ) or 0(a 2 ) , (10) 

where H = P~ + V(r) 
2m 

R 

and Hi 



2m 

pt or Hi 



?4 



m 



m 



Pi 



(11) 



for the different versions of GUP defined by equations 
(2) and (4), respectively. Thus, we see that any system 
with a well defined quantum (or even classical) Hamil- 
tonian Ho, is perturbed by Hi, defined above, near the 
Planck scale. It may be mentioned that it is often as- 
sumed that the dimcnsionless GUP parameters, ao and 
/3q are of order unity. Also, the momenta p are expected 
to be much less than M p ;c/ao or M p ;c/\//5o, because of 
the immensity of the Planck scale. Therefore the higher 
order terms which have been left out are sub-leading and 
unimportant. As for the leading order terms that have 
been retained, the resulting effects although small, may 
be perceptible, which as shall see in later sections, may 
be further enhanced by various kincmatical factors 2 . 



III. PARTICLE PRODUCTION: PARAMETRIC 
RESONANCE 

We review and closely follow the analysis of [4, 7], and 
consider an oscillating scalar field 



ip(t) = hcos(6t) 



(12) 



2 Recently, three of the current authors, without assuming the 
strengths of «o,/3o o. priori, derived bounds on the linear GUP 
parameter [25] from experimental accuracies of various quantum 
systems such as the Landau levels, the Simple Harmonic Oscil- 
lator, the Lamb Shift, and the tunneling current in a Scanning 
Tunneling Microscope (STM) etc. From the first three quantum 
phenomena, the upper bounds on ao were found to be 10 23 , 10 17 
and 10 10 , respectively. The first bound leads to a length scale 
bigger than the electroweak length scale, and so although not 
inconsistent with it, does not lead to any new verifiable scale. 
The second and third bounds could signal a new and intermedi- 
ate length scale between the electroweak and the Planck scale. 
It was also found that even for ceo & 1, one might be able to de- 
tect quantum gravitational corrections in an STM. Additionally, 
the bounds on the quadratic GUP parameter [20, 21], for the 
Lamb Shift, Landau Levels and STM were found 10 18 , 10 25 and 
10 10 , respectively. In other words, the bounds derived for the 
linear GUP parameter arc more stringent than those derived for 
the quadratic GUP parameter. It is possible that suitable astro- 
physical or cosmological observations may also put some useful 
bounds on the above parameters. 
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coupled to another scalar field x representing the matter, 
via an interaction 



r 1 2 



(13) 



Then the evolution equation for \ is nothing but the 
Mathieu equation, of the form 



x- 



+ ujI (l + h cos (2lu + e)t ) X = 



(14) 



where the argument of the cosine, 9 = 2uj +e, is so chosen 
to produce the strongest parametric resonance (PR) via 
the ft-term. We shall assume h,e/ujQ << 1, and retain 
terms only to leading order in h, e. Assuming a solution 
of the form 

X(t) = ao(t) cos (^t) + b (t) sin (^-t) (15) 

with a ~ e St)t , b ~ e St)t , substituting in Eq. (14), using 
ao ~ eao, &o ~ (thereby ignoring ao <~ e 2 ao and 6o ~ 
e 2 6o terms), identities such as cos ^4 cos B = i[cos(A + 
B)+cos(A-i?)] and cos A sin B = |[sin(A+.B)-sin(A- 
B)], and ignoring weaker resonance terms of the form 
cos(n§i) and sin(n|i) (neN> 1), we get 



a s 



Mo 



bo ( hujp 
2 V 2 
ao / hujp 
Y I" 



+ e 



— e 



= 
= 



Solving the above two equations, we obtain 




(16) 
(17) 

(18) 
(19) 



Note that < R < oo, with R = 1 corresponding to 
ao = fro and R = (oo) corresponding to &o = (ao = 0). 
Thus, when s € K, the solution given in Eq. (15) grows 
exponentially in the so-called Instability Region of huio 
given by the parameter range 



(20) 



where /(x, x) is an arbitrary non-linear function. As we 
shall see later, A is suppressed by powers of Planck mass, 
and its effects may only show up at very high energies 
and very small length scales. We will thus treat this 
term pcrturbatively. Once again, we assume Eq. (21) 
has a solution of the form 



x(t) 



a(t) cos y-t) + b(t) sin ^ -t 
where a ~ e st , b ~ e s 



(22) 
(23) 



with a = ao + Aai, 6 = 6o + A6i, s = So + Asi . 
Thus we can substitute the 'unperturbed' solution given 
by Eq. (15) in the RHS of Eq. (21), simplify again us- 
ing trigonometrical identities, and retain only the leading 
order resonance terms to write 



A/(X, X) = Aw sin ( -t ) /i(a , b , s Q ) 



+Aa; cos ( -t ) / 2 (ao,&o,s ) 



This modifies Eqs.(16-17) to 



as 



hug 



a fhuj 
bs+ 2{— 



2 /2 • 



(24) 

(25) 
(26) 



Now we make a further assumption about the smallness 
of A, namely that A <~ e 1+p , < p < 1 and retain terms 
up to C(e 2 ) to obtain 



A 



bo 



huo \ h ( huj o , 



• (27) 



Comparing with the leading order expression s 2 = Sq + 
2AsqSi we get 



2s si 



A / h^o _ \ _ h ( l^o 
bo V 2 7 a V 2 6 



(28) 
(29) 



Thus to determine whether PR is enhanced (si > 0) or 
diminished (si < 0), one would need to find expressions 
for /i , f 2 for specific models of non-linearity. 



This is the phenomenon of Parametric Resonance [4] . 



PARAMETRIC RESONANCE WITH GUP 



IV. PARAMETRIC RESONANCE WITH 
NON-LINEAR TERMS 

We introduce a generic non-linearity in the RHS of Eq. 
(14) of the following form 

X + uj 2 (l + h cos [(2w + e)t] ) X = A/(x, x) (21) 



We start with a GUP-modificd Hamiltonian in one di- 
mension of the form 

H-^v M .i + v&) + !^ { , (30) 

where k — ±1. This incorporates the two versions of 
GUP presented in Section II, with n = 4 and n = 3, 



respectively, (cf. Eq. (11)). The first of the equations of 
Hamilton, and its inverse, both to 0(a n ~ 2 ) are given by 



dH po una 11 2 „ -, 

X = jt- = — + Po 

opo m m 



(31) 
(32) 



Po = mx — nna n 2 (mx) n 1 . 
Then the second Hamilton equation p = —dH/dx gives 

m X [l- Kn(n-l)a n - 2 m n - 2 X n - 2 ] = (33) 



.. dV dV .„„ 
or, mx+-^- + Ki—x = 

dx dx 



(34) 



where k\ = Kn(n — l)a n 2 m n 2 . Thus for the time de- 
pendent harmonic oscillator, we get from Eq. (34) above 



X + wg(l + fccos (0t)); 



(35) 



where we have ignored terms 0(n\h). Thus , comparing 
Eqs. (21) and (35), we obtain Xf = - KiLu 2 xx n ~ 2 ■ Fur- 
ther, we will make the identification A = a™~ 2 . Next we 
write the solution given in Eq. (22) as 



X = ce 4 2* + c *e * 2 * 
X = Ae^ + A*e~ i i t 



(« + 



(36) 
(37) 

(38) 



from which the quantity xx" 2 that appears in the RHS 
of Eq. (35) reads 



= E( n A*("- 2 -J') (cei(2j-n+3)|t + (40) 



J=0 



c * e i(2j-n+l)|^ 



(41) 



Setting 2 j — n + 3 = ±1 and 2j — n + 1 = ±1 to extract 
the e ±l 2 4 , i.e. the dominant resonance terms, we see j = 
(n - 2)/2, (n - 4)/2 and j = n/2, (n - 2)/2, respectively. 
It is evident that n must be even and for the rest of the 
paper, we assume that this is the case 3 . Collecting these 
terms, simplifying, using a — sa, b — sb, and replacing 
{a, b, s} by {ao, bo, so} by noting that the above term is 
multiplied by the small parameter A in the equation of 



motion, and finally comparing with Eq. (24) we get 



h 

h 

where £ 





n 




n — 


2 


= ial 


n 




n — 


2 



; (4s 2 + 9 2 )R + R{As 2 - O 2 ) - 4s o 

{As 2 + e 2 ) + {As 2 -e 2 ) + AR SQ e 

K 2 (l + -R 2 ) 

16 

n- 2 



(42) 
(43) 
(44) 



and k 2 = -nn{n - l)m n - 2 io ( ) | A\ n ~ 4 



(45) 



which using Eq. (29), implies in terms of the ratio R, 

2soSl = M(^ +e ) [(rf _ 1)( _^ K+e2) 



+ (As 2 - # 2 )) - 8Rs 6 



(46) 



We use the last equation to summarize the sign of s\ 
(Note that sq is positive) 





sign of si 
[it > 0) 


sign of si 
(k < 0) 


R= 1 


+ 




R = 


+ 




R = GO 


+ 




R^l 


± 


± 



Thus, we see that the initial conditions on the matter 
field (via R) and the GUP one is considering (via the 
sign of k) determine whether there is an increase in the 
exponent of the matter field or not. The various auxiliary 
variables that were introduced (such as A, C, I etc) do 
not play any role in it. 



A. Instability Region 

Now setting s = in Eqs. (25) and (26), we see that the 
PR occurs when the modified instability region is given 

by 



--huj t— < e < 

2 60 2 ao 



(47) 



whose width is thus 

Ae = hu>o 



X 
aoR 

A(l + R 2 ) 



= huj - X 

H 

= kujQ — 2Xhujo&£a 2 ) 
= hu \l - 2X0£al 



(Rf 2 - h) 
s o 0£ 



(48) 



Higher order (weaker) resonance terms can of course arise for n 
odd. 



It follows that parametric resonance is maintained in the 
presence of GUP. However the region of instability in- 
creases when k > and vice- versa, i.e. it once again 
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depends on the GUP under consideration, as well as the 
amplitude and oscillation frequency of the matter field, 
which may contribute to further enhancement of this re- 
gion (but not on R). However, since such an increase or 
decrease would be proportional to inverse powers of the 
Planck mass, it may turn out to be too small to have an 
observable effect at present, although with more accu- 
rate experiments and improved observations, it may be 
detectable in the future. In the next section, we re-do 
the analysis for an expanding universe. 



VI. PARAMETRIC RESONANCE IN AN 
EXPANDING UNIVERSE WITH NONLINEAR 
TERMS 

We adopt the procedure outlined in [4, 6] . We first write 
the equation for the matter field in the presence of the 
inflaton field ip and in an expanding background [1], to- 
gether with the non-linear terms described by Eq. (30) 
as 

X + 3Hx + wg [1 + /icos^o + e)t] X = A/(x, *), (49) 

where H = a /a is the Hubble parameter and a(t) is the 
scale factor for FRW spacetime. The re-definition 



Y = a 3 / 2 



X 



implying 



X = a-*l\Y-°-HY) 



X = a- 3/2 (Y +^H 2 Y -3HY -^HY) (52) 

substituted into (49) yields the following equation for the 
function Y 



(50) 



(51) 



Y + uj 2 Y = XP[Y,Y,a] (53) 



where 



Lo 2 {t) = ujI(1 + h cos 6t)-^H 2 -^H 

with 9 = 2uj q + e , (54) 
and P = a 3 / 2 f[a-^ 2 Y,a-^ 2 (Y -\hY)\ . (55) 

Next we consider the following cases. 



A. Case when H and A = 

To gain some insight we first revisit the case where 
there is no nonlinearity. We assume a solution of the 
form 



in which the effects of expansion are included in the scal- 
ing function Co( a )- Substituting into Eq. (53) with its 
RHS set to zero, one obtains 

{A 1 + A 2 ) cos ^tj + (Bi + B 2 ) sin (^tj + 0(6 2 , h 2 ) = 

(57) 

with 

M = (d 6 + u 2 c + ^c - ^)Co (58) 



(59) 



B 1 = {-coO + <4do - -±-d a - ^)Co 
A 2 = 2c Co + cqCo + doCo# - jH 2 c (o - ^Hc ( (60) 
B-2 = 2d Co + d C - c Q ( 9 - ^H 2 d ( Q - ^rf Co(61) 

If Eq. (57) is to be justified, the coefficient of the sine 
and cosine must vanish. In addition, in order to ensure 
resonant behavior, we further set A\, A 2 , B\, and B 2 
separately equal to zero. Thus, up to to order 0(e,h), 
the coefficients A x and B x reduce to the Eqs. (16) and 
(17), i.e., 







(62) 



Ai = 2s d + c (-jp - 

B 1 = -2s co-rfo(^ + e) =0 (63) 

where s is the characteristic ('unperturbed') exponent 
defined in Eq. (19). On the other hand, A 2 = and 
B 2 = yield 



*»(§M§)+*(§)»-(M*)« 



'Ci 



Co 



l H2 + l H 
9 „o 3 • 







which can be rewritten as 

^.| (44) + (fc (|),_(J a . + |A)«,_ ( M) 

<ii«<»)-«>(i) e -(i fl2+ ^)* = « 

The two equations above may be combined to give (set- 
ting c = C e s °' and d = D e Sot in which C and D 
are constants) 

|(e 2sot Co)-(^ 2 + ^)Co e 2s "^0. (66) 

This equation may be solved following the procedure 
given in [6]. For a cosmic-time scale factor of the form 
a(t) oc t q , we have 4 

H + H 2 = q(q - l)t~ 2 = a{q)t- 2 , H = qt' 1 , (67) 



Y(t) = coCo(a) cos ( -t J + doCo(a) sin I -t ) (56) 



4 It is noteworthy that ct(q) > for q > 1, and that for de Sitter 
expansion, a(oo) = oo. 



6 



then Eq. (66) becomes 

Co + 2s Co - P{q) t- 2 (o = 



(68) 



where we have introduced the coefficient (i(q) = \q 2 + 
|a(g). If we use the fact that so is of order e and h, then 
we are allowed to approximate the above equation to 

Co - /?(<z)^ 2 Co * . (69) 
The solution to this equation takes the form 

Co(t) = c i aa+viw) + C2 ^(i-7iw) 

3q 2 — 3q 

= C\t~ + C 2 t— . (70) 
With this solution for (o(t), the resonant solution reads 

'0 



Y ± (t) = CoCo(t)e ±Sot 



cos ( -t I T Rs'm I -t 



(71) 



where R stands for c?o/co- 

Following the normalization chosen in [11], Co has been 
determined to be Co = \J~-^§- By comparing the solution 

given by Eq. (70) to the solution associated with the 
nonresonant case as given in [6] , we infer that the second 
exponent in Eq. (70) must be neglected. Thus, the ap- 
propriate exponent for the resonant case is the first one. 
The requirement that Co[a(i)] = 1 when a{t) = 1 im- 
plies the normalization Co = a(t) 3 / 2 . Therefore, the full 
resonant solution in an expanding background behaves 
as 



Y ± (t) 



s±« 



COS 



-t =F -Rsin I -t 
2 / \2 



with the characteristic exponent s±(t) defined by 
s±(t) 



±Sn + Tt lnt - 



(72) 



(73) 



Clearly, the stability band width will decrease after tak- 
ing the expanding background into account. In partic- 
ular, following Eq. (20) the bounds of the Instability 
Region in a flat spacetime are 



and 



while when the background is an expanding one as de- 
scribed above the bounds are modified as follows [4, 6, 9] 





= +\ 



(75) 
(76) 



Now we are in a position to generalize the above calcu- 
lation to the nonlinear model governed by Eq. (53), and 
find the analogue of Eq. (66). 



B. Case when H ± and A ^ 

We reconsider equation (53) with the nonlinear terms 
included 

Y + uj 2 Y = XP[Y, Y, a] (77) 
and thus, the solution will now be of the form 
'9 



Y(t) = c(t)C(a) cos ^-t j + d(t)C(o) sin ^-t j (78) 

with c = Co+Aci, d = do+Xdi, s = s +\s 1 , ( = Co+ACi • 
Next we substitute the unperturbed solution given in Eq. 
(56) in the right-hand side of Eq. (77), and as in the 
previous sections perform the trigonometrical approxi- 
mations and retain term of order 0(h). We thus write 



XP[Y,Y,a] = Awosinj-t) Pi [c , d , s , Co, a] + 

Xuj cos \^tj P 2 [c ,d ,s ,(o,a} -(79) 
Then we obtain up to order 0(h, A) 



(B - XuoPi) sin \-tj+(A- Xlj P 2 ) cos ( -t 



with 





+o(e 2 ,h 2 ,x 2 ) = 


(80) 




B = B 1 +B 2: A = A 1 +A 2 , 


(81) 


e modification of Eqs. (58) - (61) to read 






/ - n 7 , u 2 h , d0 2 , A 
= (-c9 + u 2 Q d - ^-d - — )C 


(82) 


A x 


= (dO + Lu 2 c+^-c-—)C 


(83) 


B 2 


= 2d( + d( - c(0 - ^H 2 d( - ^Hd( 


(84) 


A 2 


= 2c( + c( + die - ^H 2 c( - ^Hc( . 


(85) 



Following the same methodology of the previous section, 
namely demanding resonant behavior, one can set 



B x = XujaPx 

Ai = Xcv P 2 

B 2 = 

A 2 = . 



(86) 
(87) 



(89) 



From the first two equations and up to order C(e, h, A) 
it follows that 



2sc + <^ + e)=-A(g), 
e )=A(g). 



2sd + c(^-, 



(90) 
(91) 
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We may now assume that A <~ e 1+p with < p < 1, and 
keep terms up to order 0(e 2 ), we then find 

2 ^ = icb(^ + e )t RP '- ft l (92) 

which reduces, for Co = 1 when a = 1 to Eq.(29). 
Rewriting equations A 2 — and B 2 = as in the previ- 
ous section yields 

|(e^)-(^ 2 + fffW rt = (93) 



where s is defined by 



s = s + Asi 



(94) 



and si is given by Eq. (92). At this point, it should be 
noted that s± is time-dependent, so is s. Equation (93) 
may be written out as 

t+(2\s 1 t + 2 S )(- (^H 2 + 1h)c = (95) 

giving at order 0(A°) and 0(A), respectively, 



(96) 



,i +2.s„o - ( |h 2 + |h] = -2(;.s J ./ + .s- lK() .(;<)7) 

The solution to the first equation has been determined 
in the previous section to be 



Co(t) = ^ . 



(98) 



The task now is to substitute Eqs. (98) and (92) into 
Eq. (97), and ignore terms of order e as they are small 
compared to the others. We employ Eq. (92) in order to 
evaluate Si but for this purpose we also need to specify 
the form of functions Pi and P 2 . This would determine 
explicitly function £i and thus the complete resonant so- 
lutions in an expanding background where nonlinearities 
are present will be of the form 

Y ± (t) oc e [ ±S0+ ^ ln t+X ^ + W^[co^ TRsm^-tj] (99) 

where si, Co, and Ci are an known for a given nonlinearity 
expressed by Pi and P 2 . 



C. Parametric Resonance in an expanding 
background with GUP 

In this subsection we consider the nonlinearity to be the 
GUP term, namely, 



XP = - Kl cjlxX n - 2 



(100) 



with 



X = a-*' 2 Y, X = a-*l\Y - ±HY) . (101) 

In order to find the expressions of P\ and P 2 defined in 
Eq. (79), one can follow the steps outlined in section V, 
and simply make the following replacements 



-3/2 



<, b 



-3/2 



d(, 



1 



-3/2 



A -> B = a' 3 ' 2 (A - '^H{c - id)(^j 



(c-id)C, 
(102) 



with A = \[{cC+cC,+dC,^)-i{dC,+dC,-cC,^). K should be 
stressed that the quantity "a" that appears on the right- 
hand side of the above replacements is the scale factor of 
the expanding universe. Thus the expressions of Pi and 
P 2 read 



P 1 =l Q a-»t 2 cld 

— (Asl + 2 )R + R{Asl-e 2 ) 



Pi 



with 



to a 
n 

n-2 
co 

AO) _ 



9/2 3^3 

(As 2 + 9 2 ) + (As 2 - 9 2 ) + ARs 9 
to 



4 0) (i + i? 2 ) 



16 



(103) 

(104) 
(105) 



and «2 = —Kn(n — l)ra 



n-2 



LO 



n/2 



B \ n -\106) 



where B a is the quantity B as defined in Eq. (102) but 
in which c, d, ( have been replaced with cq, do, Co, respec- 
tively. 

Note that these two expressions reduce respectively to 
Eq. (42) and Eq. (43) when we take the limit H = 
0, a(t) = 1, and Co( a ) — 1- Furthermore, the expression 
of s\ as given in Eq. (29) now takes the form 



si 



£ oa - 9 / 2 c 2 Co 2 /^ 
As V 2 

+ (4s 2 - 2 )) - 8Rs 6 



(R 2 -i)(^(A4 + o 2 ) 

(107) 



where n can be equal to 4 or 3 depending on the version 
of GUP under consideration. In Section V was shown 
that when n is even one obtains the dominant resonance 
terms, while when n is odd one gets the higher order 
resonance terms which are weaker. Thus, employing Co = 
C e Sot , a(t) — t q , H = qt^ 1 , Co = , and keeping only 
terms of order O(e ), Eq. (107) for n = 4 becomes 

t 



inm 2 u)l , 4 s / huj 

s i = ; \ K ~ V hr 



4s 



e\C 2 e 2ts ° 



Accordingly, Eq. (97) is now written as 



(108) 



Ci - f3(q)t- 2 d 



(109) 
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The solution of Eq. (109) for £1 together with the expres- 
sion (108) of si will determine the form of the solution 
(99). 

Once again, from Eq. (108) one can infer about the sign 
of Si, which as in the case of a static background, depends 
only on the GUP considered (sign of k) and the initial 
conditions on the matter field (via R). 





sign of si 
(«>0) 


sign of si 
(«<0) 


R > 1 




+ 


R< 1 


+ 





It is evident that s\ = for R = 1. 



D. Instability region 

In this case, the instability region is obtained by setting 
the exponent s of the complete resonant solutions in an 
expanding background with nonlinearities (Eq. (99)) to 
zero. The exponent is 



s = ± So + %nt + \( Sl + ^-) 

It (fit 

and setting s — and using Eq. (19), one gets 



Ci 



(110) 



with the instability region given as usual by Ae = e max — 
tmin- A number of comments are in order here. First, 
it is easily seen that switching off GUP by setting A = 
reduces Eqs. (Ill) and (112) to the ones relevant for an 
expanding universe without non-linearities, cf. Eqs. (75) 
and (76). Similarly, setting A = and q = reduces 
the above to ordinary PR, and the corresponding range 
Eq. (20). Finally, although by setting q — one should 
in principle recover the width given by Eq. (48), (by 
finding the solution for (\ in Eq. (109) and using the 
expression for si in Eq. (108)), this is seen much more 
easily by setting s = in Eqs. (90) and (91), leading to 
the instability band 



huio 



A 



c R(o 



(RP2-P1) 



(113) 



which on using Eqs. (103) and (104) for Pi and P 2 , 
respectively, reduces to 



Ae = hivo — 



R 



4(l + i? 2 )s 6 



= huj a - 2\hujo6e clCoa~ 9/2 
= huj \l - 2A6»4)CoCo a ~ 9/2 • 



In this case too, the increase in the instability region 
depends on the GUP parameter, as well as parameters 
related to the expansion of the background, which may 
further magnify the GUP effect. Still the effect may re- 
main small and unobservable due to powers of inverse 
Planck mass (via A), although the situation can change 
with better experiments and observations in the future. 
It is evident that setting a(t) = 1 and Co = lj an d also 
replacing c — > a , reduces the above to Eq. (48). 



(114) 



VII. CONCLUDING REMARKS 



In this paper, we studied the effect of GUP on para- 
metric resonance in post-inflation preheating in a static 
as well as an expanding background, and showed that 
depending on the exact form of the GUP and initial 
conditions, the phenomenon of parametric resonance 
and the corresponding instability band can increase, po- 
tentially resulting in higher rates of particle production 
and an early termination of the above. We believe the 
inclusion of GUP takes into account (at least partially) 
remnant Planck scale effects on the reheating of the 
universe. It would be interesting to study the effects of 
back-reaction of the produced particles in our set-up. 
It would also be interesting to apply our approach to 
string-inspired models, for which nonlinearities enter 
the matter field equation via the Born-Infield term, viz, 
C <~ 1 — \J\ — x 2 . This would perhaps shed some light 
on the inflationary reheating theory in the context of 
extended objects, e.g, D-branes. Last, since our matter 
field Hamiltonian (Eq.(9)) is non-relativistic, it might 
be worthwhile studying its relativistic generalization 
(e.g. via a GUP modified Klein-Gordon equation), to 
have a better understanding of the massless limit among 
other things. Note that our general formalism, including 
Eqs. (21) and (49) remain well suited for such generaliza- 
tions, as well as for a large class on non-linear corrections 
to parametric resonance and particle production. We 
hope to report on these elsewhere. 
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